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ABSTRACT 



Aims. To investigate the temporal evolution of temperature, emission measure, energy loss and velocity in a C-class 
solar flare from both an observational and theoretical perspective. 

Methods. The properties of the flare were derived by following the systematic cooling of the plasma through the response 
functions of a number of instruments - the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI; >5 MK), 
GOES-12 (5-30 MK), the Transition Region and Coronal Explorer (TRACE 171 A; 1 MK) and the Coronal Diagnostic 
Spectrometer (CDS; ~0.03-8 MK). These measurements were studied in combination with simulations from the 0-D 
Enthalpy Based Thermal Evolution of Loops (EBTEL) model. 

Results. At the flare on-set, upflows of ~90 km s^^ and low level emission were observed in Fexix, consistent with pre- 
flare heating and gentle chromospheric evaporation. During the impulsive phase, upflows of ~80 km s~^ in Fexix and 
simultaneous downflows of ~20 km s^^ in He I and Ov were observed, indicating explosive chromospheric evaporation. 
The plasma was subsequently found to reach a peak temperature of ^13 MK in approximately 10 minutes. Using 
EBTEL, conduction was found to be the dominant loss mechanism during the initial ~300 s of the decay phase. It was 
also found to be responsible for driving gentle chromospheric evaporation during this period. As the temperature fell 
below ~8 MK, and for the next ~4,000 s, radiative losses were determined to dominate over conductive losses. The 
radiative loss phase was accompanied by significant downflows of <40 km s~^ in Ov. 

Conclusions. This is the first extensive study of the evolution of a canonical solar fiare using both spectroscopic and 
broad-band instruments in conjunction with a hydrodynamic model. While our results are in broad agreement with 
the standard fiare model, the simulations suggest that both conductive and non-thermal beam heating play important 
roles in heating the fiare plasma during the impulsive phase of at least this event. 

Key words. Sun: Flares - Hydrodynamics 



1. Introduction 

The temporal evolution of most solar flares can be di- 
vided into two distinct phases. During the impulsive 
phase temperatures rise to >10 MK via direct heating 
below the reconnection site in the corona and/or the 
process o f chromospheric evapora tion due to accelerated 
particles (|Kopp fc Pneumaul 119761 ) . Chromospheric evap- 
oration driven by accelerated particles can be classified 
in one of two ways - exp losive or gentle ([Fisher et al.l 
Il985t iMilligan et all l2006allbl ) . Explosive evaporation oc- 
curs when the flux of non-thermal particles impact- 
ing the chromosphere is greater than a critical value 
(~3x 10^° ergs cm~^ s""'^) and the chromosphere cannot dis- 
sipate the absorbed energy efficiently enough. The plasma 
is forced to expand into the corona as hot, upflows of hun- 
dreds of km s~^ and simultaneously into the chromosphere 
as cooler downflows of tens of km s~^. Beam driven gentle 
evaporation occurs when the non-thermal flux is less than 
~10^° ergs cm^^ s^^. Under these circumstances, the chro- 
mospheric response is efficient in radiating the absorbed 
energy. Gentle evaporation can also be driven by a down- 
ward heat flux from the corona. In both cases, plasma rises 
slowly (tens of km s"'^) upwards into the loop. Once the 



energy release has ceased, the hot plasma returns to its 
equilibrium state during the decay phase. The cooling pro- 
cess begins with thermal conduction as the dominant loss 
mechanism due to the high temperatures present. As the 
temperature decreases and the radiative loss function be- 
gins to increase, rad iative cooling becomes more efficient 
(jCulhane et al.|[l970f) . Finally, the "evaporated" material 
drains back towards the solar surface, returning the system 
to equilibrium. 

There have been a wealth of studies that focus on hy- 
drodynamic modelling of these heating and cooling mecha- 
nisms (e.R.[An tiochos & Sturrock, 1978; Fisher ct al., 1 
' Doschek et al.l. |1 983; Cargill, 1993; Klim chuk fc Cargil 
20011: iReeves fc Warren. ' ,2002; .Bradshaw fc Cargil 



Klimchuki. 120061; IWarren fc Winebargerl 1200 



ISarkar fc Walshl . l2008f) . For example, |Reald (12007D con 
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ducted an analysis of the details of s tellar flares using 
the Palermo-Harv ard theoretical model (jPeres "eTaniigsa 
iBetta et al.lll997t) . This paper fully describes the cooling 
timescales and plasma parameters of flares in terms of their 
phases, including an investigation of the thermal heating 
function. However, these, and most other theoretical 
results were not compared to observations. The majority 
of investigations that make this comparison concentrate 
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Fig. 1. Evolution of the GOES C3.0 flare observed on 2002 March 26 between 15:00 UT and 16:30 UT. RHESSI 6- 
12 and 12-25 keV images, integrated over two minutes (the time shown ±1 minute) are shown in the top two panels. 
The contours from Ov and Fexix observed at 15:09 UT are overplotted in red and yellow respectively. Fexix (8MK), 
Fexvi (2.5 MK), TRACE/171 A (1.0 MK) and Ov (0.25 MK) emission is shown in the next four panels. The dotted 
box overplotted on the Fexix images represents the region of the loop used in the temperature and emission measure 
investigation. The bottom panel shoes the GOES 1-8 A, RHESSI 3-6, 6-12 and 12-25 keV lightcurves. RHESSI was in 
eclipse until 15:15 UT and passed through the South Atlantic Anomaly between 15:35 UT and 15:58 UT. The vertical 
dotted lines (and corresponding arrows) on the GOES plot represent the start and end times of the CDS rasters above. 
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on broad-band instruments an d utilis e very simple models. 
For example ICulhane et al.l (|1994D compared Yohkoh 
observations to an over-simplified power-law cooling curve. 
lAschwanden fc Alexander! (|2001h compared broad-band 
observations to a model that considers a purely conductive 
cooling phase f ollowe d by a purely radiative cooling phase. 
IVrsnak et al.l (|2006l ) conducted a similar study, again 
concentrating on a broad-band observations and a simple, 
independent cooling mechanism model. 

The work presented in this paper aims to improve 
on previous studies by comparing high resolution obser- 
vations over a wide range of temperatures to a detailed 
theoretical model. Observations of a GOES C-class so- 
lar flare were made with several instruments, including 
the Transition Re gion and Coronal Explorer (TRACE, 
iHandv et al.l Il999( l. the Reuven Ram aty High Energ y 
Solar Spectroscopic Imager (RHESSI; iLin et all 120021) . 
GOES-12 and the C oronal Diagnostic Spectrometer (CDS; 
iHarrison et"aI]ll995^ . There are many advantages to using 
spectroscopic data in conjunction with broad-band observa- 
tions. The identification of emission lines are, for the most 
part, well documented and therefore individual lines can 
be isolated for analysis. Also, material as cool as 30,000 K 
can be observed simultaneously with emission at 8 MK. 
Furthermore, it is possible to carry out velocity, tempera- 
ture and emission measure diagnostics over a wide range 
of temperatures for the duration of the flare, significantly 
improving the scope of the analysis undertaken. These ob- 
servations were compared to a highly efficient 0-D hydro- 
dynamic model - t he Enthalpy Based T hermal Evolution 
of Loops (EBTEL: iKlimchuk et all[200l . 

The combination of this extensive data set and the new 
modelling techniques enables a comprehensive analysis of 
the heating and cooling of flare plasma to be performed. 
Section [5] describes the observations of this flare and the 
diagnostic tools used. Section [3] contains a summery of the 
theoretical models and Sect. H] lays out the results of this 
study. The conclusions and future work are discussed in 
Sect, m 

2. Observations and data analysis 

This investigation concentrates on a GOES C3.0 flare that 
occurred in active region NOAA AR9878 on 2002 March 26 
close to disk centre (-92", 297"), beginning at ~15:00 UT. 
The CDS observing study used (FLARE_AR) focused on 
five emission lines spanning a broad range of temperatures. 
The rest wavelengths and peak temperatures of each of the 
fines and those of the RHESSI, GOES and TRACE pass- 
bands are given in Table [U where the quoted temperatures 
refer to the maximum of the response function. Each raster 
consists of 45 slit positions, each 15 s long, resulting in an 
effective cadence of '^11 minutes. The slit itself is 4" x 180", 
resulting in a 180" x 180" field of view. 

Figure [T] shows the evolution of the flare in multiple 
wavelengths. The top two rows show the loop top source 
observed in RHESSI 6-12 and 12-25 keV energy bands, with 
the Ov and Fexix (15:09 UT) contours overplotted. The 
next four rows of this figure shows emission observed in 
Fexix, Fexvi, TRACE 171 A and Ov. The GOES 1-8 A 
and RHESSI 3-6, 6-12 and 12-25 keV lightcurves are shown 
at the bottom of the figure. 

As Fig. [T]shows, at ~15:00 UT, before the main impul- 
sive phase of the flare began, evidence of low level Fexix 



loop emission was observed. By ~15:09 UT, the footpoints 
were seen in O v while the Fe xix loop top emission contin- 
ued to brighten. At 15:16:40 UT, when RHESSI emerged 
from eclipse, a thermal looptop source was observed in both 
6-12 and 12-25 keV energy bands, with a corresponding 
temperature of '--^15 MK. By 15:21 UT the dominant emis- 
sion had cooled to ~8 MK. At this time, a bright "knot" can 
be seen at the top of the loop. Such features have been ob- 
served in the past and hav e not been readily explained (e.g. 
iDoschek fc Warrenl . l2005f ). During the early decay phase 
(~15:32 UT), the loop is seen to cool into the Fexvi tem- 
perature band (~2.5 MK) and by ~15:44 UT, the plasma 
has cooled to <1 MK, as seen by TRACE and in Ov. 

Although RHESSI was in eclipse for the majority of the 
impulsive phase of the flare, the observed continued rise of 
the 6-12 keV lightcurve after emergence from night, im- 
plies that the peak of the soft X-rays (SXR) was observed. 
However, while a hard X-ray (HXR) component was ob- 
served at this time, it is not believed to be the peak of 
non-thermal emission. 

Table 1. Rest wavelengths and temperature of emission 
lines and bandpasses used in this study. 



Ion 


Ao [A] 


Temperature [MK] 


He I 


548.45 


0.03 


Ov 


629.80 


0.25 


Mgx 


625.00 


1.2 


Fexvi 


360.89 


2.5 


Fexix 


592.30 


8.0 


TRACE 


171 


1.0 


Instrument 


Range 


Temperature [MKJ 


GOES 


0.5-4A & 1-8A 


5-30 


RHESSI 


3 keV-17 MeV 


^5 



2.1. Temperature and emission measure 

The temperature and emission measure evolution for 
this flare was determined by analysing lightcurves from 
RHESSI, GOES, CDS and TRACE. 

The RHESSI spectrum, shown in Fig. [31 was analysed 
over one minute at the peak of the 6 -12 keV energy band . 
Following previous studies (e.g. ISaint-Hilaire fc Bend . 
.20021 ), the data were fitted with an isothermal model at low 
energies and a thick-target model up to ~30 keV. The thick 
target component yielded a low energy cutoff of 17 keV and 
a power law index of 8.2. The isothermal fit to lower en- 
ergies determined a temperature and emission measure of 
^--^13 MK and ^1x10** cm~^ respectively. A non-thermal 
electron flux of ^7x10^ ergs cm~^ s""'^ was also calculated, 
approximating the footpoint area from He i and O v obser- 
vations. Since it is probable that the HXR peak occurred 
before this time, this value is taken to be a lowe r limit to the 
maximum non-thermal electron flux. FoUowing' White et al.l 
(2005), the filter ratio of the two GOES passbands produced 
a temperature profile, giving a peak temperature of 10 MK 
and an emission measure of 4x10'*^ cm~^. 

CDS and TRACE 171 A observations were integrated 
over a fixed area described by the bright "knot" , seen in 
Fexix at 15:21 UT in Fig. [TJ This area is highhghted by 
the dotted box in the Fe Xix images. The lightcurves of this 
region were analysed for each raster and TRACE image 
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function for that particular temperature and density. The 
Fig. 2. Linear plot of the hghtcurves of emission observed intensity, /, was calculated using: 
by GOES, CDS and TRACE. The data points were fit using 

a spline interpolation. The vertical represent the peaks of /" m\ 2 iTr ^1 mNrnTir ^/ m\ 2 i-Tr /i\ 

the hghtcurve fits. ^ ^ / G{n^^T)nldV = G{n^,T)EM = G{n^,T)nljV (1) 



during the course of the flare (Fig. ^ . The loop is believed 
to consist of multiple magnetic strands. It is assumed that 
the majority of strands within the region of this "knot" are 
heated almost simultaneously at the time of the HXR burst. 
A small number of strands can be heated before or after 
this time, producing a multi-thermal plasma. However, this 
small region is approximately isothermal at any one time. 
For the remainder of this paper, any reference to the "loop 
apex" refers to the area defined by the dotted box on the 
Fexix images in Fig.[TJ The Hei line was not used for the 
determination of temperature as it is optically thick and 
neither He i nor TRACE were used in the emission measure 
analysis. 

To obtain the thermal evolution of this flare, the time 
of the peak of each lightcurve was assigned the associated 
temperature mentioned in Table [T] The uncertainty in the 
temperature measurement was taken to be the width of the 
appropriate contribution function. The motivation for this 
assumption is indicated in Fig. ID Figure [4^ shows the con- 
tribution functi£n_forMgx calculated using CHIAN TI ver- 
sion 5.2 (jPere et al.. Ill997t iLandi fc Phillipsll2006[ ). Figure 
[4)d shows a theoretical temperature evolution for a flare 
with a maximum of ^20 MK and Fig. U): shows the cor- 
responding density evolution. The lightcurve (Fig. HJi) was 
calculated by taking the density and temperature at each 
time step, and calculating the value of the contribution 



for constant volume, V ^ temperature T, filling factor, /, 
electron density, and emission measure EM. G{n^,,T) 
is the contribution function calculated for constant den- 
sity, nJ3 and variable temperature, T. During the impul- 
sive phase, as the corona is heated, the temperature jumps 
by ~2 orders of magnitude, passing quickly through the 
peaks of the contribution function and density curves. As 
this happens, a sharp peak is seen in the lightcurve at 
'^lOOO s. As the plasma cools, a broader peak is seen in 
the lightcurve, relating to emission from the looptop. This 
occurs at ~2600 seconds (dashed lines. Fig. [Hd, c, d), a 
time corresponding to a temperature of 1.25 MK, the tem- 
perature at which the contribution function is maximised 
(dotted lines. Fig. 3^, b). Therefore, the lightcurve peak 
seen in the cooling phase of the flare can be attributed to 
passing through the maximum of the contribution function; 
i.e., at 2600 s, the plasma is emitting at 1.25 MK. 

Under the assumption of an isothermal plasma, the 
emission measure was obtained at the times of the 
lightcurve peaks for Ov, Mgx, Fexvi and Fexix, follow- 
ing Eqn. [T] We have not included TRACE in the analysis 
of EM due to the ill-defined instrument response of the 



^ A study of the contribution function dependence on density 
was carried out. It was found that varying n, by two orders 
of magnitude resulted in little or no change to the contribu- 
tion function. This density was assigned a constant of value of 
10^^ cm-^ 
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Fig. 4. Panel (a) shows the Mgx contribution function cal- 
culated from CHIANTI. The peak is highlighted with a dot- 
ted vertical line. This occurs at 1.25 MK. Panel (b) shows 
the simulated temperature profile calculated by the EBTEL 
model (see Sect. |3] for details). The dotted horizontal line 
denotes the temperature of the contribution function peak 
from panel (a) , and the dashed vertical line is the time this 
temperature is reached. Panel (c) shows the simulated den- 
sity profile, and panel (d) shows the simulated lightcurve 
for Mgx, calculated from Eqn.[TJ 



TRACE 171 A band. The uncertainty in calculating EM 
using RHESSI was found to be approximately 50% using a 
thermal fit to the spectrum. This is in agreem ent with the 
value s found by the RHESSI instrument team (jMcTiernanl 
[2006h . 

For CDS, there are a number of factors to consider. 
These in include uncertainties in the intensity of the line, 
the contribution function and the CDS calibration. While 
the uncertainty in measuring the line intensity is small 
for str ong lines such as those used in our study, typically 
~10% (|Del Zanna et al.ll200lD . a consideration of the con- 
tribution functions FWHM results in an uncertainly in the 
EM of approximately 30%. The CDS calibration is also 
known to be good to within 15-20% (jBrekke et'al.|[2000l) . 
The GOES instruments is known to have limited ability 
for making accurated measure ments of both temperature 
and EM. Accepted values from lGarcial ()2000D . give the EM 
uncertainty to be 10%. Considering these factors, the com- 
bined photometric error for all instruments was taken to be 
50%. 



Fig. 5. The line profiles of the Fexix emission line for the 
right footpoint during the impulsive phase (a) and the de- 
cay phase (b). Note the increasing intensity of the Fexii 
blend at ~592.8 A (jPel Zanna fc Masod[2^5h during the 
decay phase. The solid vertical line represents the rest wave- 
length while the dashed vertical line is the centroid of the 
line. 



2.2. Velocity 

The Doppler shifts at both footpoints were calculated for 
the duration of the flare using the five CDS emission 
lines, with uncer tainties of ±10 km (jBrekke et al.lllQQTl : 
iGallagher et al.llTQQa) . The rest wavelengths for Hei, Ov, 
Mg X and Fe xvi were calculated from a region of pre-flare 
quiet Sun. For the Fexix line, a detailed analysis of the be- 
haviour of wavelength as a function of time was conducted 
to establish the rest wavelength. Centroid wavelengths from 
both before the SXR rise and from late in the flare were 
averaged and corrected for both heliocentric angle and an 
average inclination of 44° to obtain the rest wavelength. 
Fig. [5] shows the Fe xix line profiles for the right footpoint 
during the impulsive (a) and decay (b) phases. 



3. Modelling 

3.1. The Cargill model 

Follo wing lAntiochos fc Sturrockl ()1976[ ). ICargilll (|1993l 
Il994l ) presented a model that considered a fiare that is cool- 
ing purely by conduction for a time Tc, followed by purely 
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radiative cooling for a time t^. The cooling time-scales were 
given by: 



4 X 10 



-10 



y5/2 ' 



and 



3fcT 



(2) 



(3) 



where A(T) is the radiative cooling function (see 
iRosner et al.l Il978l for further details) , L is the loop half 
length, and all other variables have their usual meaning. 

The time and temperature at which the cooling mech- 
anism dominance changes, r* and respectively, can be 
calculated as follows: 



7/12 



and 



= T,; ( - 

Tr 



-1/6 



(4) 



(5) 



where T* = T{t = t*) and subscript denotes initial 
values. These parameters were calculated for the flare under 
investigation and compared to the results of the 0-D hydro 
model. 



3.2. The EBTEL model 

Enthalpy Based Thermal Evolution of Loops (EBTEL) 
model is a OD model that simulates the evolution of the 
average temperature, density , and pressure along a single 
strand (jKlimchuk et al.|[2008D . calculating a single value of 
each of these quantities at any given time. This is a reason- 
able representation since temperature, density, and pressure 
are approximately uniform along the magnetic field, with 
the exception of the steep gradients in the transition region 
at the base of the loop. 

As its name implies, EBTEL takes explicit account of 
the important role of enthalpy in the energetics of evolving 
loops. Under static equilibrium conditions, less than half of 
the energy deposited in the corona is radiated directly. The 
rest is thermally conducted down to the transition region, 
where it is radiated away. Under evolving conditions, chro- 
mospheric evaporation occurs when the transition region 
cannot accommodate the downward flux, or, if the flux is 
insufficient to power the transition region radiative losses, 
condensation occurs. 

EBTEL equates an enthalpy flux with the excess or 
deficit heat flux. Kinetic energy is ignored because the flows 
are generally subsonic, except perhaps in the earliest times 
of an impulsive event. Another assumption made is that 
the radiative losses from the transition region and corona 
maintain a fixed proportion at all times. EBTEL has been 
compared with sophisticated ID hydrodynamic models and 
found to give similar results, despite using 4 orders of mag- 
nitude less computing time. 

EBTEL allows for any temporal profiles of both direct 
plasma heating and non-thermal particle acceleration. The 
effects of the non-thermal electron beam are treated in a 



highly simplified manner. It is assumed that all of the en- 
ergy goes into evaporating plasma. This is reasonable for 
gentle evaporation (Fisher ct al. 1985), but for explosive 
evaporation, some of the beam energy will go into a plug 
of downflowing and radiating plasma deep in the chromo- 
sphere. Thus the actual energy of the beam, as inferred 
from RHESSI observations for example, is greater than the 
beam energy used in the EBTEL s imulation. For a comple te 
description of this model, refer to iKlimchuk et al] (|2008t ). 

The flare loop is almost certainly composed of many 
strands that are heated at different times. However, the ob- 
servations suggest that most of the strands are heated in ap- 
proximately the same way and at approximately the same 
time (i.e. during the HXR burst), so the flare was modelled 
as a single monolithic loop. Nonetheless, some strands are 
expecte d to be heated both bef ore and after this main bun- 
dle (e.g.. lKlimchuk et al.l . l2006h . and this will result in some 
deviations between the model and observations. 

The pre-flare conditions included a temperature of 
0.3 MK, an initial density of 5x10^ cm""^ and an emis- 
sion measure of 4x10'*'^ cm^'^. Input values were, where 
possible, constrained by observations. The loop length was 
determined from magnetic field extrapolations of the re- 
gion (Conlon, priv. comm.). The filling factor could not 
be constrained by data due to the poor spatial resolution 
of the high temperature instruments. Since it was assumed 
that the heating function is associated with the HXR burst, 
the majority of which was not observed, the shape of the 
heating function was inferred from previous observations of 
HXR bursts and the slow rise of the GOES SXR lightcurve. 
The most appropriate heating function was Gaussian in 
shape. The non-thermal electron flux was constrained by 
that found by the lower limit calculated from RHESSI ob- 
servations and the wid th was i nferred from the deriv ative of 
the SXR flux (jNeuper t 1968; Dennis fc Zarrolll993t) . While 
the direct heating rate was not constrained by observations, 
it was assumed to have the same width as the non-thermal 
heating flux and to occur at the same time. Due to the sen- 
sitivity of the model parameters to cooling timescales, the 
cooler data points (e.g. Fexvi, Mgx, TRACE and Ov) 
were critical in constraining the parameters. The ranges 
of acceptable parameter values are shown in Table [21 The 
values obtained from observations are shown, along with 
the parameter values used in producing the results in Sect, 
m The range of parameter values shown correspond to the 
maximum and minimum values that produce an acceptable 
fit to data. The ratio of the heating components (i.e. direct 
to non-thermal) is also shown for the best fit parameters 
presented in Sect.Hl where the equivalent direct energy flux 
is given by the volumetric heating rate divided by the half 
loop length. 



4. Results 

Combining the observations from the different instruments 
used for this study, and the results from EBTEL, the heat- 
ing and cooling phases of this flare can be comprehensively 
described. These results are presented in Figs. [6] and [71 
These figures show the evolution of the flare through the de- 
pendence of temperature, emission measure, energy losses 
and velocity, as discussed in Sects. [2TT] and [2?2l 
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(A) (B) (C) (D) 




15:00 15:15 15:30 15:45 16:00 16:15 16:30 

Start Time (26-Mar-02 14:47:02) 



Fig. 6. The EBTEL temperature evolution (solid line) that best reproduced the observations of Sect. I2.1l is shown in the 
first panel, along with the heating function whose parameters are described in Table [2] (dashed line, arbitrarily scaled to 
show position and width). The second panel shows the corresponding model and observed emission measure evolution, 
with the data points corresponding in time to those in the first panel. The third panel shows the conductive and radiative 
losses throughout the flare on solid and dashed lines respectively. The fourth panel shows the GOES 1 ~ 8 A lightcurve. 
The bottom panel shows the velocities of the right footpoint for the five CDS lines and that of the left footpoint in He i 
and Fexix. The velocities were very similar at both footpoint in all other lines and so were omitted for clarity. The 
EBTEL simulated velocity for Mgx is represented by the thick black line. The dotted vertical lines correspond to the 
flare phases (A)-(D) explained in Sect. 14.21 
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Table 2. Input parameters used for modelling simulation. The parameters were constrained by data when possible and 
the ranges of parameters investigated are shown. 



Parameter Observed EBTEL 

Loop half-length [cm] JxW (3 ± 0.2) x 10'' 
Non-thermal flux 

- Amplitude [ergs cm"^ s"^] 7x10" 5 x 10**^ 

- Width [sec] ~100 100 ± 50 

- Total [ergs cm"^] ~1.7xl0^^ 2.5 x 10^°=^^ 
Direct heating rate 

- Amplitude [ergs cm~^ s~^] - 0.7 ± 0.3 

- Width [sec] - 100 ± 50 

- Background [ergs cm~"^ s~^] - < 1 x 10~^ 

- Total [ergs cm"^] - 175 ± 150 
Direct/non-thermal heating (best fit parameters) ~ 4 



4.1. Comparison of model to data 

The top two panels of Fig. [6] describe the evolution of the 
flare temperature and emission measure, respectively. The 
data points for each instrument were obtained using the 
analysis described in Sect. 12.11 The input parameters for 
EBTEL were approximated by observations and allowed 
to vary slightly until a good correlation with the cooling 
phase data was obtained. The conductive and radiative loss 
curves generated by EBTEL for the flare are shown in the 
third panel of Fig. [S) Consistent with previous observations, 
conduction was found to dominate initially, with radiation 
becoming prevalent for the remainder of the decay phase. 
Both Cargill and EBTEL found conduction to dominate 
for the first 200 - 400 s of the decay phase, with radiation 
dominating for the remaining ~4000 s, referring to Tc and 
Tr respectively. The time at which Tc ~ (i.e. the domi- 
nant loss mechanism switches from conduction to radiation) 
is ~15:24 UT in both cases. However, the temperature, 
T, = T(r*), at which this occurs was found to be ~12 MK 
and ~8 MK according to Cargill and EBTEL respectively. 
This discrepancy is due to the different approach to the 
modelling of the early decay phase. EBTEL simultaneously 
calculates the conductive and radiative losses throughout 
a flare while Cargill assumes cooling exclusively by either 
conduction or radiation at any one time. The fourth panel 
of this figure shows the GOES 1-8 A lightcurve for con- 
text. The last panel shows the velocities at the loop foot- 
points, calculated following the analysis in Sect. 12.21 The 
flow velocity at both left and right footpoints are shown for 
the coolest and hottest lines - He i and Fe XIX respectively, 
while for clarity, only the right footpoints for the remain- 
ing three lines were shown. The MgX Doppler shift simu- 
lated by EBTEL is represented by the thick black line. The 
simulations are in reasonable agreement with the observa- 
tions. Upflows are of course predicted during the evapora- 
tion phase and downflows are predicted during the draining 
phase. However, the magnitudes are generally larger than 
those observed for reasons that we do not fully understand. 
Uncertainties in the velocity zero point adopted for the ob- 
servations may account for the downflow discrepancy. 

The redshifts observed in the cooler lines during phase 
B, if real, are likely an indication of a downflowing chromo- 
spheric plug that accompanies explosive evaporation. The 
blueshifts seen in the hotter lines around 15:30 UT suggest 
that some loop strands were impulsively heated after the 
primary flare energy release. These blueshifts are expected 
to be smaller than the actual upflows because the spatially 
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Fig. 7. This shows the dependence of emission measure on 
temperature for both model and data. The different phases 
of a solar flare are marked (A)-(D). Over-plotted are the 
emission measure data-points calculated in Sect. 12.11 as a 
fimction of their temperature. 



unresolved line profile re presents a mixture of upflowin g 
and downflowing strands ([Patsourakos fc Klimchu]jl2006l ). 
The small peaks seen at later times in the simulation ve- 
locity curve are a result of the piecewise continuous form 
used for the radiative loss function. 

Figure [7| shows the evolution of the flare through the 
interdependence of emission measure on temperature. The 
data points obtained during the analysis described in Sect. 
12.11 were computed at the same instant in time for any one 
emission line or bandpass. This figure shows the heating of 
the plasma (A) followed by evaporation of hot plasma (B), 
cooling (C) and draining (D). The plasma initially cooled 
to a temperature below the pre-fiare value and asymptoti- 
cally returned to the equilibrium state due to the low level, 
constant background heating. 

4.2. Flare Phases 

Figures [6] and [7] have sections labelled (A)-(D) which refer 
to the different phases of the flare, from pre-flare heating 
to the late decay phase and are described in detail in this 
section. 
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A) 14:45-15:10; Pre-flare phase: For the majority of 
this phase, the EBTEL parameters remained at quiet Sun 
values, as Fig. [6] shows. At 15:07 UT the EBTEL temper- 
ature and emission measure began to rise. Figure [7] shows 
the steep temperature gradient and the initial gradual rise 
in emission measure. However, as the fourth panel in Fig. 
[6] shows, the GOES soft X-rays began to rise slowly before 
this. At ~15:00 UT, a small amount of Fexix emission was 
seen in the loop (Fig. [T]) and velocities of ~90 km s~^ ob- 
served in Fe xix can be seen in the bottom panel of Fig. [6] 
while all of the cooler lines remain at rest. This is evidence 
that pre-flare heating is driving gentle chromospheric evap- 
oration in a small number of strands heated before the HXR 
burst. 

B) 15:10-15:17; Impulsive phase: During the impul- 
sive phase of a flare, the standard model predicts the 
propagation of non-thermal electrons to the chromosphere 
where they heat t he ambient plasma , caus ing it to rise 
and fill the loop (jKopp fc PneumanI [l975)- Upflows of 
^80 km s^^ in Fexix and simultaneous cool, downflows 
of ~20 km s^^ in Hei and Ov were observed and shown 
in the bottom panel of Fig. [6l A non-thermal electron 
flux of ^7x10^ ergs cm^^ s^^ was determined between 
15:16:30 and 15:17:30 UT. This is slightly lower than the 
3x10^*^ ergs cm~^ s~^ required to drive explosiv e chro- 
mosph eric evaporation (jFisher et al.l [1985; Millig an et al.l 
l2006al) . As such, this value is taken to be a lower limit and 
that the maximum value of non-thermal flux occurred be- 
fore 15:16 UT. 

C) 15:17-15:24; Soft X-ray peak: The top panel of Fig. 
[6] shows the temperature has peaked and begun to fall and 
that the emission measure and SXRs were at a maximum 
in this phase. As the third panel of Fig. [6] shows, during 
this phase, conduction was efficiently removing heat from 
the corona and transferring it to the ch romosphere, driving 
slow upflows of hot emission (jZarro fc Le men 1988.) . These 
flows ~20 km s^^ in Fexix respectively can be seen in the 
bottom panel of Fig. [51 

D) 15:25-16:30; Decay phase: This phase is dominated 
by radiative cooling, as seen in the third panel of Fig. [B] 
Velocities in Fexix were returning to quiet Sun values. The 
modest blueshift observed at 15:32 UT in Fexix suggests 
that the line profile contains components from evaporat- 
ing strands that were heated after the main loop bundle. 
Between approximately 15:45 and 16:20 UT Mgx, Ov and 
He I showed downflows of up to ^40 km s~ ^. This implies 
loop draining was occurring (|Brosiud 120031 ). By the end of 
the simulation, all of the parameters had returned to quiet 
Sun values. 

5. Conclusions and Discussion 

This paper compares a flare observed with CDS, TRACE, 
GOES and RHESSI to the 0-D hydrodynamic model 
EBTEL. Early in the impulsive phase of the flare, evi- 
dence of 8 MK emission and '^90 km s~^ upflows sug- 
gest pre-flare gentle chromospheric evaporation. During 
the impulsive phase, hot upflowing plasma at velocities of 
~80 km s""'^ and cool downflows of ~20 km s^^ imply ex- 
plosive chromospheric evaporation. Around the time of the 
soft X-ray peak, conduction was found to be highly effi- 
cient. Upflowing plasma at velocities of ~20 km s^^ are ob- 
served in Fexix, suggesting conduction driven gentle chro- 
mospheric evaporation. The cooling timescales modelled by 



ICargil]l(|1993lll994f ) were tested against EBTEL and proved 
to agree reasonably well. The dominant cooling mecha- 
nism was found to switch from conduction to radiation at 
^--^15:24 UT in both models. However, the temperature at 
which this occurs {~8 MK and ~12 MK for EBTEL and 
Cargill respectively) does not agree. This is as a result of 
the simultaneous cooling by conduction and radiation for 
EBTEL versus the exclusive cooling phases predicted by 
ICargilll (|1993| . Il994[ ). The late decay phase of the flare is 
dominated by radiative cooling. Downflowing plasma ob- 
served in He I, Ov and Mgx during the late decay phase 
provides evidence of loop draining. 

By tracking the behaviour of this flare as it cooled 
through the response functions of the many instruments 
and emission lines, the evolution of the temperature and 
emission measure could be assessed. This evolution was 
then recreated using the EBTEL model, providing precise 
details, such as the cooling timescales and mechanisms, 
that cannot be easily obtained from data. For this par- 
ticular flare, since the HXR burst was not fully observed, 
the details of the heating function could be also estimated 
from simulations. The description of the flare using both 
data and model allows for a much greater understanding 
of flare dynamics. The nature of these explosive events re- 
main somewhat ambiguous, however further studies of this 
nature, will help to improve the understanding of them. 

The ratio of the heating functions were investigated. It 
was found that the observations were best modelled when 
the plasma was heated approximately equally by direct and 
non-thermal mechanisms. This implies that both of these 
processes are vital during the flaring process and that flares 
may not be energised primarily by non-thermal particles, 
as previously believed (lBrown| | 1971l). This i s in agreement 
with recent results found bv lMiUiganl (gOOl). There it was 
shown that a non-thermal electron beam is not necessar- 
ily required to obtain the high-temperature, high-density 
material we see in flares. However, it should be noted that 
the EBTEL value of the flux of non-thermal electrons re- 
quired for equal heating is below the critical value for ex - 
plosive evaporation hypothesised by iFisher et al.l (jl985l ). 
This can be explained by the over-simplifed treatment of 
non-thermal particles by EBTEL. This requires caution 
for a flare of this nature, where it is evident that non- 
thermal particles play an important role. The model as- 
sumes that all non-thermal energy is used for evaporating 
plasma upward into the loop. However, this may not be 
entirely true. It is well known that a very small fraction of 
this energy is used to produce bremsstrahlung radiation (1 
part in 10^). It may also be possible that a more signifi- 
cant amount is used to force plasma down int o the chromo- 
sphere and to power chr omospheric emission (jWoods et al.l 
l2004lAllred et al.ll2005l ). However, despite the approxima- 
tions made by EBTEL, such as the homogenous nature of 
the loop or the disregard for the location of energy deposi- 
tion, the temperature and emission measure curves repro- 
duce observations very well. It is computationally efficient, 
running a complete simulation in a matter of seconds. 

This paper has established a method that will be ap- 
plied to the analysis of future events. For this case, the data 
was manually compared to theoretical model. However, the 
fitting of the parameters together with model comparison 
techniques is currently being investigated using a Bayesian 
technique for simulating values from the posteri or distribu- 
tions of the parameters (jAdamakis et al.ll2008[ ). The pur- 
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pose of this analysis is to statistically optimise the model 
parameters within boundaries set by observations. This ap- 
proach will be used when comparing theoretical models to 
future data sets. The authors intend to carry out an inves- 
tigation of flare hydrodynamics using the improved cadence 
and extensive spectral range of the Extreme ultraviolet 
Imaging Spectrometer (EIS) on board Hinode. Combining 
these data with RHESSI spectral fits will vastly improve 
observations and allow for even more accurate modelling. 
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